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Abstract 

Simulating the evolution of the Human Immunodeficiency Virus (HIV) epi- 
demic requires a detailed description of the population network, especially 
for small populations in which individuals can be represented in detail and 
accuracy. In this paper, we introduce the concept of a Complex Agent Net- 
work(CAN) to model the HIV epidemics by combining agent-based modelling 
and complex networks, in which agents represent individuals that have sexual 
interactions. The apphcabihty of CANs is demonstrated by constructing and 
executing a detailed HIV epidemic model for men who have sex with men 
(MSM) in Amsterdam, including a distinction between steady and casual 
relationships. We focus on MSM contacts because they play an important 
role in HIV epidemics and have been tracked in Amsterdam for a long time. 
Our experiments show good correspondence between the historical data of 
the Amsterdam cohort and the simulation results. 

Key words: Multi- Agent Systems(MAS), Complex Nctwork(CN), HIV 
epidemics, scale-free. Men who have Sex with Men(MSM) 



1. Introduction 

Understanding the underlying dynamics in Human Immunodeficiency 
Virus (HIV) epidemics is a crucial pubhc health issue, unfortunately how- 
ever addressing specific problems in small populations is difficult because 
individuals need to be modelled with detailed social behavior. Traditional 
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mathematical methods greatly simplify both the disease dynamics and the 
population networks, however extending them to more detailed models is 
intractable. 

In particular, whether such methods also apply to networks of small size, 
and thus to many real- world biological or epidemiological applications, is 
still an open question[l]. Many networks of relevance to epidemiology may 
be of relatively small size, among which social contact networks on sexually 
transmitted diseases for a small group (e.g. men who have sex with men, 
namely MSM, within a city or town) are representative. 

Modelling the HIV epidemic is difficult because the true incidence^ of 
the HIV/AIDS-epidemic is uncertain since many people may be unaware of 
their infection. Secondly, HIV progression has a very long asymptomatic 
period which makes studies of the actual infection spreading a very compli- 
cated task [2] . Finally, the various routes of infection and the inhomogeneity 
of the involved population pose additional challenges to understanding the 
underlying knowledge of HIV epidemics. 

Multi-agent systems (MAS) and complex networks (CN) are often used 
separately to model and simulate epidemics; however, whether existing mod- 
els, which typically focus on large populations, can address epidemics among 
small groups (~ 10^-10^) are still questionable and seemingly need further 
vahdation. In terms of agent-based modelling, Teweldemedhin developed 
an agent-based bottom-up modelling approach for estimating and predict- 
ing thcX spread of the HIV in a given population [3]; Xuan developed an 
extended Cellular Automata simulation model to study the dynamical be- 
havior of HIV/AIDS transmission by incorporating heterogeneity into agents' 
behavior [4]. In terms of complex network modelling, Bai discussed a sexual 
network spreading model for HIV epidemics [5]; Sloot proposed a new way 
to model HIV infection spreading through the use of dynamic complex net- 
works, with the time evolution of the network vertices modelled by a Markov 
process [2]. 

In this paper we present the CAN approach for simulating epidemics in 
small networks to great detail. The CAN is a hybrid approach in which multi- 
agent systems and complex networks are the basic methods of modelling 



^Incidence is a measure of the risk of developing some new condition within a specified 
period of time. In this paper we use the incidence rate to denote the number of new 
cases per unit of person-time at risk. 
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epidemics on individual and population scales respectively. Using the CAN 

approach, we simulate a relatively detailed model of the HIV endemic among 
MSM in Amsterdam. This model includes a distinction between steady and 
casual relationships, which is regarded as an important aspect by others [e.g. 
6] . We compare the results to the Amsterdam Cohort Study (ACS) historical 
data. 

This paper is organized as follows. In Section 2 we review MASs and 
complex networks and introduce the CAN approach. We describe the HIV 
epidemic model in detail in Section 3 and simulation implementation in Sec- 
tion 4. In Section 5 we discuss a case study and its simulation results. 

2. Complex Agent Network 

The CAN approach takes the MAS and CN as basic methods for mod- 
elling and simulating epidemics on individual and population scales. Agents 
contain specified personal information on an individual scale, whilst complex 

networks emphasize the relationship dynamics among the agents on a popu- 
lation scale. We review both concepts and the combination of them in this 
section. 

2.1. Multi- agent Systems 

A MAS is a system composed of multiple interacting intelligent agents, 
and can manifest self-organization and complex behavior even when the 
individual strategies of all the agents are simple. Heylighen defined self- 
organization as "the spontaneous emergence of global coherence out of the 
local interactions between initially independent components" [7]. 

The agents in a multi-agent system have several important characteristics 
[8]: Autonomy: the agents are at least partially autonomous. Local views: 
no agent has a full global view of the system, or the system is too complex for 
an agent to make practical use of such knowledge. Decentralization: there is 
no controlling agent (or the system is effectively reduced to a reductionistic 
system) [9]. 

Applied to HIV epidemics, a MAS represent a specified human commu- 
nity and viruses propagate along social contacts. Each individual has his 
own progression of infection and simple rules to choose partners. Each agent 
only knows the information about himself and his partners instead of epi- 
demiological statistics. There is no dominating agent in the community to 
control the spread of viruses. 
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2.2. Complex Networks 

The study of complex networks is inspired largely by the empirical study 
of real-world networks such as computer networks and social networks. A 
network is a set of items, which we call vertices, with connections between 
them, called edges [10]. In the context of network theory, a complex network 
is a network with non-trivial topological features that mostly do not occur 
in simple networks such as lattices. 

The complex network theory aims to do three things [10]. Firstly, it aims 
to find and highlight statistical properties, such as path lengths and de- 
gree distributions, that characterize the structure and behavior of networked 
systems, and to suggest appropriate ways to measure these properties. Sec- 
ondly, it aims to create models of networks that can help us to understand 
the meaning of these properties, how they came to be as they are, and how 
they interact with one another. Finally, it aims to predict what the behavior 
of networked systems will be, on the basis of measured structural properties 
and the local rules governing individual vertices. 

For HIV epidemics, an at-risk population can be depicted as a network 
with vertices representing persons and edges representing social contacts. 
Over these relationships, viruses arc propagated from one individual to an- 
other and thus give rise to the epidemic prevalence. Findings in sociology 
or epidemiology suggest that social contact networks are scale-free networks 
[11]. Therefore we create networks with vertices following power-law degree 
distributions, and then simulate HIV propagation by evolving networks tem- 
porally and spatially. 

2.3. Complex Agent Networks 

The benefits of combining a MAS and a CN is twofold. 

Firstly, the edge forming and rewiring mechanism of a complex network 
with a non-trivial topology can regulate the interactions among agents in a 
MAS. In this way we can keep these interactions from being either purely 
regular or purely random as we usually see in a MAS. Furthermore, the 
recent development of statistical methods for quantifying complex networks 
is a way to analyze the infiuence of network topologies on epidemics [e.g. 12]. 

Secondly, local agent information can be combined to infer global statis- 
tics which can be fed back into the simulation dynamically, resembling real- 
world epidemics more closely. For example, if the HIV prevalence becomes 
very high, individuals may be less prone to unsafe sexual interactions with 
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Infer global statistics 




Figure 1 : Complex Agent Networks 

other individuals whose serostatus they do not know.^ This interaction of 
both the individual and population levels is depicted in Figure 1. 

To combine the MAS and CN we map agents and their interactions in a 
MAS to vertices and edges in a CN, respectively, as shown in Fig. 1. Agents 
provide personal heterogeneity and local infection progression on an individ- 
ual level; networks benefit from continuously distilling values of epidemic- 
related parameters from individuals and then mimic the virus propagation 
on a population level. 

2.4- Application of CAN to HIV epidemics 

Fig. 2 shows an application reference procedure of applying the CAN 
approach to HIV epidemics, which consists of four steps, i.e. data collection, 
modelling, simulation and result manipulation. 

The CAN-based HIV epidemic studies commence with data collection 
from e.g. Clinical Databases, Internet Data Mining, Literature and Ques- 
tionnaires. 

Fed with the data, models are constructed by utilizing the features of 
both the MAS and CN, which are depicted by two composite blocks. The 
"My4S"' composite block deals with societal individuals and their infection 
details. An Agent stands for an individual and supports specialized inheri- 
tance. As shown in Fig. 2, an inherited block MSM Agent manages MSM' 



^Here we assume that individuals somehow find out such global statistics, e.g. through 
media. 
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Figure 2: Application reference procedure of Complex Agent Network 



infection properties and status transitions by using HIV Progression and 
Treatment blocks. A MSM Agent also support inferring values of parameters 
for a ^^CN" . The Treatment block considers factors related to HIV ART or 
HAART'^ (e.g. drug types, patients' adherence and drug resistances [13]). 
The ^^CN" composite block manipulates network initializing, reshuffling and 
measuring through the Network and Network Measurement blocks. In ad- 
dition, the ^^CN" provides network rendering support for the Visualization 
block of simulation. 

The Visualization and Checkpoint blocks can be invoked by the Simu- 
lator according to specified purposes, to enhance interactive exploration of 
simulations. 

The Statistics and Analysis blocks are applied to simulation results and 
then enforce the GUI Update of network rendering and parameter plotting. 

3. Modelling HIV Epidemics 

3.1. Progression of HIV infection 

The progression of HIV infection from seropositive to AIDS can broadly 
be divided into three stages: primary infection (PI), asymptomatic period 
(AP) and AIDS. Although there are much variations amongst patients, on 



■^ARV, Anti-Retroviral Therapy; HAART, Highly Active Anti-Retroviral Therapy, is 
the combination of at least two different classes of antiretroviral drugs. 
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average, transmissibility is the highest during the primary infection which is 
bcUcvcd to be associated with high plasma HIV RNA levels and continued 
risky sexual activities [14, 15]. 

Correspondingly each agent representing an at-risk individual, possesses 
a disease manager to conduct its HIV infection progression temporally. The 
manager conducts an irreversible "negative — ^ PI — ^ AP — ^ AIDS" order. 
Note that, we exclude the case that an individual with AIDS can with the 
right treatment go back to the asymptomatic phase, because in general pa- 
tients are most likely diagnosed and treated before reaching the AIDS stage. 

3.2. Construction of contact networks 

HIV is transmitted mainly through sexual contacts, intravenous drug uses 
(IDU), mother-to-child transmissions and transfusion of contaminated blood 
products. In this paper we focus on MSM contacts which play an important 
role in HIV epidemics and have been tracked in Amsterdam for a long time 
[e.g. 16, 17]. We model these contacts as an undirected, scale-free complex 
network. 

The sexual contact network starts with Nq MSM, each of which is assumed 
to be sexually active with an age between 15 and 65 based on the age limits 
of the ACS samples [17]. The network is assumed to be scale-free [11, 18], 
whilst the degree k of each vertex follows a degree distribution ~ k~'^ to 
represent his yearly number of partners. Let us choose 



where 7 and po are constants and C is a normalizing factor. 

After obtaining the degree sequence of vertices, edges (partnerships) are 
formed by using the so-called configuration model [10] . A degree sequence is 
a set of n values of the degrees ki of vertices i = 1, n, from the given power- 
law distribution. According to the configuration model, the probability that 
a vertex is chosen as an end of a random edge is proportional to kpk- 

Different from most previous studies on epidemics in static networks, 
our networks evolve with simulation time steps, which is driven by the HIV 
epidemic dynamics depicted in Section 3.3 and 3.4. Further implementation 
will be discussed in Section 4.2. 




for A; > 1 
for A; = 



(1) 



7 



3.3. HIV propagation in networks 

The HIV propagation in networks can become considerably difficult when 
more details are taken into account. We provide a formula of calculating the 
transmission probability, considering risk behavior, treatments, (steady or 
casual) partnerships and transmissibility at different infection stages. 

Similar to the commonly used Susceptible-Infected-Removed (SIR) epi- 
demiological model, our model sets individuals in three discrete states, i.e. 
susceptible, infected or removed (caused by deaths). The infected state com- 
prises primary infection, asymptomatic period and AIDS stages (see Section 
3.1) which are related with different statistical transmission probabihty {TP) 
per action. These TP per action are represented by TPpi, TP^p and TPaids 
for PI, AP and AIDS stages, respectively. Individuals at stage AIDS are 
assumed to have no contacts, so that TPaids = 0. 

Each individual Vi has ki partners, among whom at most one steady 
partner is allowed. This steady partnership lasts for years based on statistical 
data. With the presence of a steady partner, an individual is likely to reduce 
his risk behavior with all other casual partners, and this possible reduction is 
represented as a partnership factor. Two steady partners have many actions 
per year, while two casual partners have one action per year. 

For a vertex V^, the statistical transmission probability per action is de- 
noted by TPi. Then this transmission probability per action is adjusted 
according to the stage an individual is currently at, and the treatment re- 
duction factor Ft which is derived from treatment effectiveness. 

For an edge Eij as a partnership connecting and Vj, the risk behavior 
factor is denoted by Fr.. , which is related to the partners' attitude towards 
risk; the number of actions per year over this partnership is denoted by N^ij 

In addition, the partnership factor denoted by Fp. is subjected to whether 
the individual Vj has a steady partner, and it will influence all incident casual 
partnerships. 

Then the transmission probability from Vi to Vj per action is given based 
on parameters aforementioned 

PTA,,^Fp.*FR^.*{FT,*TPi) (2) 

Because the PI stage lasts around 3 months (0.25 year) and falls short 
of one time step of a year, we divide TPpi into TPpi i and TPpi 2 as two 

transmissibility at PI stage. Hence if Vi is at the PI stage, we also divide 
PrAij into PrAij,! and PTAij,2, and Na^^ into A^Ai^,i and iVAi,-,2- Then 
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Pta,,, 2 = Fp^ * Fr^^ * {Ft, * 

The transmission probability from Vi to Vj per year is 



r 



HIV negative Vi 




PI Vi and casual Eij 
PI Vi and steady ii^ij 
AP Vi and casual E'jj 
AP Vi and steady E^jj 
AIDS Vi 



(3) 



The 



probability of Vi infected by positive partners per year is 



(4) 



m=l 



where ki is the number of partners of Vi. 

Note that TPj, Fjv and Fp. are varying with the status of Vi and time, 
and so are Fr.. and A^^.^. with the status of Eij and time. Therefore each 
susceptible vertex gets infected with a probability given by Equation 4 at 
each time step, which contributes to the dynamics of HIV propagation. 

3.4- Demographic influences 

In this section, we discuss the demographic principles that regulate the 
conservation or deletion (substitution) of individuals and relationships. 

The vertices in networks are initialized randomly with an age from a range 
of 15 to 65, according to the age limits of the ACS samples [17]. Each edge is 
initialized with an expected duration according to the statistical real-world 
partnership duration. At each time step, an individual gets an increment of 
age and a partnership gets an increment of duration, mimicking the reality. 

If one of the following three conditions applies, then an individual will be 
removed from a network: 

(1) The individual's age exceeds 65, the upper age Umit of the ACS sam- 
ples. 

(2) The infection has developed to the AIDS stage (the individual will be 
removed at the neoct time step). 
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(3) In case of geographic migration (a fraction of randomly selected indi- 
viduals will be removed from a network. This fraction is left as tunable, e.g. 
1% assumed in our model.). 

Susceptible individuals will be added to the network as substitutes for the 
removed individuals at each time step, to keep the population size constant. 
All these newly added individuals will join the epidemic contact network 
according to the principles given in Section 3.2. 

An edge will be kept for the next time step if: 

(1) It is a steady partnership and has not expired. 

(2) A previous casual partnership is likely to be conserved with a proba- 
bility (e.g. 0.2 assumed in our model). 

4. Simulation Implementation 

We implement specific functionalities of the HIV epidemics simulation, 
such as an advanced scheduling mechanism and interactive visualization, 
using the JUNG[19] and MASON[20] hbraries. 

4.I. Functionality 

The specific functionalities of the CAN simulation contains the following 

components: 

(1) Model interpretation 

The model described in Section 3 is interpreted into an executable simu- 
lation supporting the population and individual scales. 

On a population scale, the formation and reshuffle of the edges are driven 
by the dynamics of the networks. An edge is formed based on the degrees 
and local status of two potential partners, as well as the (steady or casual) 
type and duration of their partnership. 

On an individual scale, heterogenized agents are described with different 
HIV progression stages (Susceptible, PI, AP and AIDS) which are influenced 
per time step by the disease progression, risk behavior and treatment. 

Furthermore, the interpreted executable simulation consists of modules 
which we call operators (see Section 4.2). We embed these operators into 
either agents or networks, and then schedule them using a simulator. 

(2) Simulator 

We implement a flexible scheduling mechanism, interactive visualization 

and a checkpoint mechanism into the simulator. The scheduling mechanism is 
double-precision discrete-event-and-priority driven and supports complicated 
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Table 1: Simulation operators 



Operator 


Scale 


Priority 


Description 


Infection 
Operator 


Population 





Infect susceptible individuals 
and propagate HIV. 


Local 

Progression 

(3l)erator 


Individual 


1 


Advance HIV disease progression 

of individuals. 


Demographic- 






Enforce demographic and 


Reshuffling 


Population 


2 


reshuffling influences on network 


Operator 






evolutions. 








Statistics approaches and 


Statistics and 


Population 




visualization updates can be 


Visualization 


3 


added and configured according 
to pre-specified outputting rules. 



user interactions with the simulator. This gives support to interactive explo- 
ration of the running simulations and, for instance, allows checkpoint/restart 
of long running simulations. 
(3) Statistics and analysis 

We can take specific epidemic-related statistics and analyze their results, 
besides general statistical properties e.g. path lengths. For example, we take 
incidence per person-year as a measure of HIV epidemics in Section 5, by 
dividing the number of newly infected persons every year by the population 
size. 

4-2. Operators 

The simulation operators with their annotated priorities arc shown in 
Table 1. These operators are invoked at each time step to evolve a network 
N[t] = {V[t], E[t]) where t is the time step. 

The Infection Operator basically deals with virus propagation on a popu- 
lation scale. For simplification, we disseminate its calculation to individuals 
by using Equations 2-4 in Section 3.3. Algorithm 1 gives the details. 

The Local Progression Operator and Demographic - Reshuffling Operator 
are presented in Algorithms 2 and 3 respectively. 
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Algorithm 1 Infection Operator 
Input: Network N = {V,E) 

Output: The evolved network N with some vertices newly infected. 

Infect susceptible vertices, 
for Vi^V do 

if Vi is susceptible then 

Judge whether Vi has any steady partnership 

Calculate PrVi using Equations 2-4 

if a generated random number < PTYi then 

Set the infection stage of Vi to Primary Infection, 
end if 
end if 
end for 



Algorithm 2 Local Progression Operator 
Input: Network N = {V,E) 
Output: The evolved network N. 

Make vertices and edges progress, 
for Vi^V do 
Vi.Kg^++ 

if Vi is infected then 

if Vi.CurrentStage.LastingPeriod < 
Vi.CurrentStagc.ExpcctcdDuration then 

Vi.CurrentStage.LastingPeriod++ 
else 

Vi progresses to next stage i.e. PI— >^AP or AP— *^AIDS. 
end if 
end if 
end for 

for Ei^E do 

if Ei is a steady partnership then 

if £^j.LastingPeriod< fii.ExpectedDuration then 

Ei .LastingPeriod-l— I- 
end if 
end if 
end for 
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Algorithm 3 Demographic-Reshuffling Operator 
Input: Network N = {V,E) 
Output: The evolved network N. 

Evolve the network demographically. 
for Vi^V do 

Remove Vi according to the rules in Section 3.4. 
end for 

Add new susceptible vertices to substitute removed ones. 

for Ei&E do 

Remove Ei according to the rules in Section 3.4. 
end for 

Reshuffle all edges according to the rules in Section 3.2. 



4.3. Simulation scheduling 

Wc design and implement a flexible scheduling mechanism with three 
main advantages. 

Firstly, a double-precision discrete-event based scheduler is available. 
With this we can schedule special events, such as the Gay Parade Day^ 
parties for MSM happening at fractional time steps. 

Secondly, all agents and operators can be registered with predefined pri- 
orities. Vertices (agents) and edges (partnerships) are scheduled firstly, (de- 
mographical and reshuffling) dynamics of networks secondly, and statistics 
and GUI updates finally. In this way, the infection and local progression 
drives the dynamic evolution of networks and the collection of statistics. 

Finally, the simulator can schedule, for special purposes, vertices and 
edges at predefined intervals. For example, we can collect incidence statistics 
every 2 years instead of 1 year. 

4.4- Visualization 

The visualization of contact networks can present an overview of the 
studied population, and explain details about a selected vertex or edge. The 
visualization serves the following purposes: 



''The Gay Parade Day is the national celebration day on the first weekend in August 
in the Netherlands. 
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Figure 3: Initialization(a) and After reshuffling(b) in the Amsterdam MSM simula- 
tion (The green, red, blue and black circles represent the susceptible group, the primary 
infection group, the asymptomatic group and the AIDS group, respectively; Solid and 
dotted lines represent steady partnerships and casual partnerships, respectively; small, 
medium and large circles represent vertices with a degree less than 50, between 50 and 
100, and greater than 100, respectively.) 



(1) To represent individuals in a population and (steady or casual) rela- 
tionships. 

(2) To arrange and display vertices and edges, in order to study the 
sorting, grouping and correlation of vertices and edges. 

(3) To present detailed properties of selected vertices or edges. 

Fig. 3 shows the visualization of the network after initialization and 
reshuffling. 

5. Case Study: the HIV epidemic in the Amsterdam MSM cohort 

We take the HIV epidemic among the Amsterdam MSM cohort 
study. 

5.1. Amsterdam cohort study 

The Amsterdam Cohort Study (ACS) of HIV infection and AIDS among 
MSM was initiated in 1984. According to the report [16], 2299 MSM have 




(a) Initialization 



(b) After reshuffling 
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been involved in the ACS until 2006. Of the 2299 MSM, 571 were HIV- 
positive at study entry and 192 seroconverted during follow-up. Yearly HIV 
inicidence for MSM over calendar years is reported in [17] to indicate the 
HIV prevalence trend among MSM. 

5.2. Parameters and their values 
5.2.1. On a population scale 

Societal-related parameters for individuals and network-related initializa- 
tion are discussed. 

Schneeberger concluded in [11] that MSM population follows a power-law 
degree distribution with a value of 7 in the interval between 1.5 and 2. We 
assume 7 = 1.6 and a maximum degree kmax = 200. For initialization, the 
network size is Nq = 2299, and the fraction of vertices with a degree of is set 
to be po = 0.01 because we assume a small portion of people in a population 
not having any contact. 

Steady relationships among MSM have a higher contribution to HIV in- 
cidence in the era of HAART because more unprotected anal sex is likely 
to take place between steady couples other that casual ones, according to 
[15, 6, 21, 22, 23]. Thus the fraction of steady relationships in a network is 
regarded as an important parameter. Based on questionnaire results given 
by [6] , approximately 50% of the ACS participants reported having a steady 
partner. So we estimate that a partnership is initialized as a steady one with 
a probability given by 

Psteady = N*{0.5/2)/{e) (5) 

where 

k 

""max 

k=0 

(e) = N*{k)/2 

N is the size of a network, (e) is the expected total number of edges in 
a network and {k) is the expected degree of vertices. According to Equation 
1 with given values of 7 and kmax, we get Pgteady ~ 0.054. Allowing for 
steady partnerships being kept until expired, the averaged fraction of steady 
relationships from our simulation results fluctuates between 0.047 and 0.055. 
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The duration of steady partnerships among Amsterdam MSM is reported 
to have an expected value of 1.5 years [6]. Thus we assume the duration of 
steady partnerships following a discrete uniform distribution W{(1,2). 

Unlike the PI or AIDS stages with fixed durations, the duration of the 
AP stage can vary with the therapy effecitiveness. So we draw the duration 
of AP from a Poisson distribution with a mean value of 13 years for a failed 
treatment, and a mean value of 22 years for a successful treatment [6]. 

The transmission probability for sexual actions other than UAI^ is as- 
sumed to be negligible [24], also, the transmission probability for an UAI 
between serodiscordant^ men is depending on whether the role of the HIV- 
negative man is receptive (URAI) or insertive (UIAI). We assume each part- 
ner of a couple takes randomly a receptive or insertive role in an action. The 
frequency of either URAI or UIAI between steady partners is 15 per year [6], 
for simplicity we assume that the frequency of sexual actions of MSM per 
year follows a Poisson distribution of V{30). At the PI stage, the frequency 
within the first 3 months and the last 9 months is assumed to follow 7^(8) 
and V{22), respectively. 

5.2.2. On an individual scale 

Epidemic parameters and the involved reduction factors are discussed. 

Of the infected individuals at the asymptomatic stage, 42% know they 
are HIV positive[6]. Since the introduction of HAART in 1996, the fraction 
of individuals that has received any form of HAART therapy after diagnosis, 
is 81% in 2006 [16]. 70% of the individuals who received therapy will be 
successfully treated[6]. In other words, individuals get diagnosed with a 
probability of 0.42, and obtain a treatment with a probability of 0.81, and 
then get successfully treated with a probability of 0.7. 

The value of TPj (see Section 3.3) is related to the infection stage of 
an individual Vi, given by TPpi^i, TPpi_2, TPap or TPaids- According to 



^UAI, unprotected anal intercourse; URAI, unprotected receptive anal intercourse; 
UIAI, unprotected insertive anal intercourse 

^Serodiscordant (sero-discordant) is a term used to describe a couple in which one 
partner is HIV positive and the other is HIV negative (Smith, Raymond. Couples. Re- 
trieved on 08-26, 2006.). Sere concordant is the term used to describe a couple in which 
both partners are of the same HIV status (ie both are HIV positive or both are HIV 
negative). 
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Table 2: The yearly relative values of risk behavior factor for HIV negative and positive 
MSM. * - The value in 1987 is an average of its adjacent years (1986 and 1988), so is that 
in 1996. ** - The baseline with value 1.00. 





1985-1986 


1987 


1988-1991 


1992-1995 


1996 


1997-1999 


2000- 


negative 


3.50 


2.50* 


1.50 


0.80 


0.90 


1.00** 


1.30 


positive 


2.80 


1.61 


0.42 


0.88 


0.78 


0.70 


1.30 



[25, 26, 27], the probability of the transmission per URAI/UIAI act at stage 
PI and AP is 0.22/0.044 and 0.011/0.0022 respectively. Therefore, TPpu = 
0.22/0.044 and TPpi,2 = TPap = 0.011/0.0022 with respect to the receptive 
or insertive role. 

ART is assumed to reduce the transmissibility by 50-90% to a moderate 
extent in [28], so that in our simulations the treatment factor Ft^ is set to 
follow a continuous imiform distribution CW(0. 1,0.5). If an individual has a 
steady partner, the partnership factor Fp^ is 0.84[see 6, Sec. 2], otherwise 1 
which means no reduction. 

The general use of ART leading to changes in risk behavior has sparked 
considerable concerns, for instance, Baggaley discussed in [15] that increases 
in risk behavior could result from increased optimism about HIV therapy 
due to the availability of HAART. In our simulations, the definition and 
estimated values of risk behavior factors FR^. is based on [29, 30], and the 
relative values are given in Table 2. 

5.3. Results and discussion 

The simulation results illustrated in Figure 4 shows the HIV incidences 
(a) and AIDS diagnosis cases (b) over calendar years. 

We used mainly the incidence for comparison and validation because of 
its relative dimension. Our simulated yearly incidences are consistent with 
the historical ACS data [see 17, Table 5], based on the conduct of a null 
hypothesis test where the simulation result of incidence does not differ sig- 
nificantly from the ACS data. Although the historical data (the dashed line 
in Figure 4a), got based on incomplete statistics (due to the involvement of 
subgroups of MSM related to substudies of ACS in years), was fluctuating 
tremendously, the hypothesis is accepted by using Chi-Square tests with a 
significance level of 0.05 (data not shown). 

As a supplement, we compared the simulated AIDS diagnosis cases with 
the historical statistics in the Netherlands in the aspect of trend. Our result 
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Figure 4: Simulation results (averaged on 24 executions) versus statistical historical data 
of HIV incidence (ACS) (a) and AIDS diagnoses cases (in Netherlands, not only in Ams- 
terdam) (b) over calendar years 

HIV Incidence Person Years (ACS Data Versus Simulation Result) AIDS Diagnoses (Netllerlands Data Versus Amsterdam Simulation Result) 




1990 2000 2010 2020 1985 1990 1995 2000 2005 



(a) Incidence (b) AIDS Diagnoses 

shows, even that it is based on a relatively small sample space with 2299 
Amsterdam MSM, a similar trend to the statistical data in the Netherlands 
(including Amsterdam as one city) over calendar years [see 31, Table B.18]. 
Both the two curves in Figure 4b approximately show a continuous increase 
in AIDS diagnosis cases until 1996 and an acute decline in 1996-1998, and 
then keeps stable from 1999 on. In this study, our simulation suggests that 
the acute decline is most likely caused by the introduction of HAART in 
1996. In the future, we could scale up the simulation of AIDS diagnosis 
cases from Amsterdam to the Netherlands, if and when the total numbers of 
MSM with different behavior features in different cities of Netherlands are 
available. 

6. Conclusions 

We proposed a flexible hybrid Complex Agent Network approach and 
validated experimentally its applicability by taking the dynamics of HIV 
propagation among Amsterdam MSM, a representative small population, as a 
case study. Our experiments showed good correspondence between the model 
results and the historical data of ACS, in addition, the results are robust 
with respect to small changes in the input parameters (data not shown). 
Therefore, the model can be adopted to predict the future trend of HIV 
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prevalence among MSM in Amsterdam, whilst allowing for easy tuning the 
values of highly concerned parameters, such as the drug effectivity, diagnosis 
proportion and risk behavior factor. 
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